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1 Introduction 

Recent years have seen a surge of activity in the area of the low x physics. Primarily 
this has been motivated by the new HERA data which has greatly extended the 
available kinematic region for Deeply Inelastic Scattering (DIS). The data now exists 
at Bjorken x as low as 10~^ and it was hoped initially that at such low x and relatively 
large one would see clearly the new class of perturbative phenomena, those that 
go under the general name of "semihard physics". First, one was hoping to see the 
unambiguous signs of the perturbative BFKL pomeron ||^ , which predicts steeply rising 
gluon distribution (and consequently DIS cross section) as a function of x at fixed Q^: 
g{x) (X (l/a;)^^-"-/''. Second, it was expected that in this kinematical region the gluon 
densities will be large enough so that the semihard shadowing effects due to gluon 
recombination will become sizable. 

In actual fact the situation turned out to be not quite so clear cut. The DIS cross section 
does indeed rise quite steeply with 1/x. It can be fit by a power of although appar- 
ently not as large a power as predicted by the BFKL formula j^]. The GLR parameter 
k which is the physical parameter for the onset of shadowing was estimated and was 
indeed found to be close to 1 within the HERA regime ||] , which would suggest observ- 
able shadowing corrections. Nevertheless, surprisingly enough all the data is described 
very well by a simple straightforward DGLAP linear evolution The rise at low x is 
then a consequence of the standard perturbative evolution of "flat" partonic distribu- 
tions from a low initial scale Qg to Q^. The shadowing corrections to this perturbative 
evolution seem to be practically absent. In this sense the higher twist corrections seem 
to be irrelevant in the HERA regime. 

Still it is widely believed that these nonlinear effects must make their presence felt when 
the partonic densities are high enough. Even if it does not happen in ep collisions at 
HERA they have a very good chance of being observed in the DIS experiments on nuclei 
at low X if and when this program takes off at HERA and even more so in the heavy 
ion collision experiments at LHC. The curious situation with the present HERA data 
only adds motivation for studying the onset of the shadowing physics. 

The shadowing regime can be approached in two ways. Decreasing at fixed (small) x 
leads to the initial growth of the cross section as Q~^. This growth gradually slows down 
and eventually (almost) stops when the unitarity bound is reached. We will refer to this 
slow down and the stop of the growth as the shadowing and the saturation respectively. 
Alternatively one can decrease a; at a fixed value of Q^. Due to the growth of gluonic 
distributions the cross section first should exhibit fast growth (powerwise according to 
the BFKL prediction), which again should slow down and saturate. Our perspective in 
this paper will be the second one, so that we will be dealing with the evolution of the 
unintegrated gluon density with 1/x. 

The physics of the nonlinear effects in DIS is basically the physics of dense partonic 
systems. This statement perhaps needs some clarification. The physical picture of shad- 
owing depends in large measure on the Lorentz frame used to describe the DIS process. 
In the rest frame of the nucleon the onset of shadowing corrections is due to the multi- 
ple scattering processes of the hadronic component of the photon (a quark - antiquark 
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pair) on the nucleon. This can be described by the extension of Glauber multiple scat- 
tering formalism to the context of QCD Q . For multiple scattering to become effective 
the partonic system does not have to be particularly dense. In this frame the onset of 
the nonlinearities (multiple scattering) is rather determined by the cross section of the 
scattering of the hadronic fluctuation on a parton in the nucleon Q. The authors of 
Ref. Q analyzed the corrections to the Glauber formula and concluded that once the 
rescattering becomes important one also has to take into account corrections due to 
rescattering of gluons in the cascade produced by the quark - antiquark pair, into which 
the photon initially fluctuates. The saturation happens in the regime where the gluon 
density in the photon is high. 

In the infinite momentum frame (IMF), where the nucleon carries all the energy prior 
to collision, the picture is different. Low x here means that one probes low longitudinal 
momentum - large wavelength fluctuations in the nucleon wave function. These long 
wavelength gluons are emitted from the valence quarks and gluons in the nucleon. In 
the standard linear evolution equations (DGLAP or BFKL) the interaction between 
these wee gluons is disregarded. However when their wave length is large enough, the 
gluons emitted by different valence partons overlap in space and interact. The first 
nontrivial effect of this interaction is the recombination process, which slows down the 
growth of the gluon density and thereby leads to shadowing 0|. The shadowing and 
saturation in this frame are both clearly effects of large partonic density^. 

Although the qualitative picture of saturation and unitarization based on the GLR type 
recombination effects is very appealing, reliable theoretical tools of dealing quantita- 
tively with finite density partonic systems are yet to be developed. The original GLR 
equation ^, |^ truncates the series in the expansion in powers of density at the first 
nonlinear term. As all such truncations it has intrinsically a very limited range of valid- 
ity, since in general one expects that when the first nonlinear term becomes important 
the higher order terms will be comparable to it in magnitude. The saturation of par- 
tonic distributions and restoration of unitarity in high energy (density) limit of QCD is 
an outstanding problem which remains unsolved although several approaches are being 
explored in the hterature l^-|l0t]. 

In this paper we continue to develop a theoretical approach to finite density partonic 
systems at low x. The main goal of this program is to derive the evolution equation for 
the gluon density at small x without assuming that the density is in any sense small. In 
the previous papers we have described the main framework of our approach and 

have discussed several aspects of this evolution. In this work we complete the derivation 
of the full nonlinear evolution equation. 

The present approach is inspired by an idea of McLerran and Venugopalan first 
formulated in the context of ultrarelativistic heavy ion collisions. The observation in 
is that there is a regime of high density and weak coupling in which semiclassical methods 

^Wc rnention here that the same conclusion emerges from the analysis carried out in the recent 
paper [o. It was shown there, using the explicit BFKL expressions for the gluon density, that for 
collision of two hadrons the shadowing first appears at low density in the frame where the two hadrons 
share the energy equally before the collision and at high density in the analog of IMF, where one hadron 
carries all the energy. The saturation again is the high density effect in both frames. 
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should apply. It was therefore suggested that the leading small x glue structure of 
the nucleus is due to the classical gluon field which is created by the random color 
charges of energetic on-shell partons. The nonlinearities of the Yang - Mills equations 
exhibit themselves already on this classical level and it is therefore possible that they 
provide the necessary saturation mechanism at low x. This approach assumes that 
the interaction of the fluctuations of the gluon field is weak. In this sense it is a weak 
coupling expansion and its validity is therefore restricted to the perturbative "semihard" 
shadowing effects. It is however nonperturbative in a different sense. In the standard 
perturbation theory, the charge density is small and the standard perturbative expansion 
is simultaneously expansion in powers of the charge density. In this language, higher 
powers of charge density appear as higher twist contributions (although there is no one 
-to - one correspondence between the twist expansion and expansion in powers of the 
charge density). The McLerran - Venugopalan (MV) method does not assume expansion 
in powers of color charge density and corresponds to resummation of a particular type 
of higher twist terms. In fact the interesting saturation effects are expected when the 
density is order aj^. This formulation is therefore naturally suited for discussion of the 
type of problems we are interested in. 

Based on this idea the approach was developed which combines the concept of the effec- 
tive Lagrangian for the low x DIS with the Wilson renormalization group resummation 
of leading log corrections to the MV approximation . The main effect of the renor- 

malization group procedure is the change in the color charge density distribution in the 
effective Lagrangian with 1/x. The RG equation that governs the evolution of this distri- 
bution is the subject of the present paper. It was shown in that in the limit of small 
color charge densities this equation reduces to the celebrated BFKL equation. In we 
have derived the general form of this evolution equation at finite color charge density. In 
the present work we calculate the "coefficients" in this renormalization group equation, 
which are in fact rather "coefficient functions" , thereby providing the last ingredient in 
the derivation of the full nonlinear evolution equation valid to leading log approximation 
at finite color charge density. We should stress, that the calculations presented here are 
only valid to leading order in as ■ The scale of as is therefore left undetermined in this 
framework and the strong interaction coupling constant is treated as a momentum inde- 
pendent constant, just like in the standard BFKL equation. Higher order perturbative 
calculations of the type of are needed to determine the appropriate scale. 

This paper is organized as follows. In Section ^ we motivate and describe the form 
of the effective Lagrangian for the physics of low x gluons in DIS. In Section || we 
describe in some detail the classical approximation to this effective Lagrangian. It turns 
out that the proper treatment of this Lagrangian requires careful specification of the 
complete gauge fixing condition, and this is also done in Section |[ We then discuss the 
first quantum corrections to the classical approximation which obviate the need for a 
renormalization group resummation, and the physical interpretation of the change of 
the color charge density distribution with the RG flow. In Section ^ we describe in detail 
the Wilson renormalization group procedure as applied to our effective Lagrangian and 
derive the general form of the RG equation. Section ^ is the central section of this paper. 
It contains the calculation of the coefficient functions that appear in the RG equation. 
Finally, Section is devoted to a discussion of our results. Several Appendices contain 



2 The EfFective Action for Low x DIS. 



5 



technical details of the calculation. 

We have tried to make this paper self contained, and therefore have included some of 
the material already contained in the earlier work ||ll|-|l3| . 

2 The EfFective Action for Low x DIS. 

Throughout this paper we will work in the infinite momentum frame, where the hadron 
moves in the positive z direction with the velocity close to velocity of light and almost 
infinite longitudinal momentum P+ — > oo. Also, we will be working in the light cone 
gauge = 0. 

Our task now is to understand the structure of the effective Lagrangian for low x DIS. 
First, it is well known that the most important degrees of freedom at low x are gluons. 
In the framework of standard linear evolution equations, the evolution of gluons in 
the leading approximation is independent of quarks, and the evolution of quarks is 
entirely driven by the gluonic distribution. We will therefore retain only gluons as our 
dynamical degrees of freedom and disregard quarks entirely. Importantly, the gluons that 
we treat as dynamical degrees of freedom are only those which have low longitudinal 
momentum, lower than some cutoff — xP'^ . Our effective Lagrangian therefore has 
to be understood as having a built in longitudinal cutoff. 

So, what is the Lagrangian that governs the interactions of the low x gluons? First of 
all, of course it must contain the standard Yang - Mills interaction term 

-Jd'x\G' (2.1) 

where G'^'^ is the gluon field strength tensor 

Cr = d'^A: - d-'A^, + gfabcA^.A^ (2.2) 

The gluons with the low longitudinal momentum also interact with the rest of the 
partons in the hadron, which have larger longitudinal momentum. We will refer to 
those partons as "fast" for notational convenience, we may think of valence partons as 
their initial representatives. This interaction certainly can not be neglected, but in the 
kinematics of IMF and in the light cone gauge it is very simple. The leading interaction 
is the eikonal vertex Jf^st; where J^^^ is the color charge density due to the fast 
partons. 

The dependence of J^^^ on x^ and x^ is very simple. First, since the wavelength of 
the fast fields is much shorter than that of the dynamical soft gluons the charge they 
produce is effectively concentrated a,t x~ — 0. Intuitively this can be understood in the 
following way. In the rest frame the valence partons are concentrated within the nucleon 
radius from the center of the nucleon. When boosted to the infinite momentum frame 
due to Lorentz contraction they are squeezed into a very thin pancake. This picture 
is a little too naive for our fast partons since some of them have much larger wave 
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length than the nuclcon radius. However as a basic physical picture it is still correct. 

We therefore have J+ oc S{x^). 

Second, we can understand the dependence by considering the (light cone) time 
scales characteristic of the problem. The relevant time scale for the low x phenomena is 
the inverse of the on shell frequency of the soft gluons: k~ oc l/k~^. The frequency of the 
fast modes is much lower p~ rx l/p'^ since their longitudinal momentum is higher, so 
that k~ /p~ oc 1/x. Therefore, as far as the soft glue is concerned the color charge source 
due to fast partons is effectively static. We are therefore led to consider the interaction 
of the type 

Si^t = A-p{xi_)6ix-) (2.3) 

The fast partons are represented in our effective Lagrangian by the surface charge den- 
sity p{x±). A hadron of course is not described by a fixed single configuration of the 
color charge density p{x±). However, the crucial point is that the structure of the fast 
component of the hadron is determined on a much longer time scale than the time scale 
relevant for the low x physics. It is fixed by the hadronic wave function, bremsstrahlung 
processes that involve fast partons, etc. Therefore, as far as the soft glue is concerned, 
there is no interference between the different configurations of p{x±). In the low x 
effective Lagrangian the hadron thus appears as an ensemble in which different configu- 
rations of p{x±) enter with some statistical weight exp{— The partition function 
for calculation of the soft glue characteristics of a hadron must therefore have the form 

J D[p, A] exp{-f [p] -'-Jd''xG^ + iSintlA, p]} (2.4) 

At this point we do not specify the form of the functional F[p]. In fact, as we shall see 
later this fmictional depends on the longitudinal cutoff A"*" which is imposed on the soft 
fields. In other words, as one considers regions of lower and lower x, F changes. The 
flow of F with the cutoff A+ a a; is described by a renormalization group equation of 
the form 

F[p] = asA[p] (2.5) 



dlnl/x 



This RG equation is precisely the evolution equation for the charge density correlators 
(and consequently for the soft glue observables) which we undertake to derive in this 
paper. 

Of course, in order to make quantitative statements about the x dependence of F we have 
to specify the initial condition for the evolution. This can be done in the perturbative 
region at not too small a value of x, where F can still be expanded in powers of p. The 
initial form of F can then be taken as 

F = J dx±dy±p{x±)n~'^{x±,y±)p{y±) (2.6) 

with 

fi{x±,y±) = S{b^) J ^ e^'=-(---^-)</.(fc^,xo) (2.7) 
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Here b i = 



is the impact parameter, S{b±) is a nucleon thickness function, xo is 



the value of x from which we start evolving F according to the RG equation and 0(fc_L , x) 
is the unintegratcd gluon density^. The relation between the parameter in the Gaussian, 
and 4>, Eq.(2.7) will become clear when in Section ^ we consider the perturbative 
calculation of the gluon structure function based on the effective Lagrangian in Section ^. 

One important question that we have not touched upon so far, is the question of gauge 
invariance of our effective action. Although we have partially fixed the gauge by the 
light cone gauge condition = 0, the action should still preserve the residual gauge 
symmetry. This residual gauge symmetry group is comprised of gauge transformations 
with gau ge fu nctions which do not depend of x~ . The naive "Abelian" eikonal interaction 
term Eq.(2.3) does not preserve this gauge symmetry. The relevant generalization takes 
the form 



5*11 



1 



(fxj_dx S{x )p°(xj^)trTaW_oo,oo[^ ]{x ,xj_) 



(2. 



Here Ta are the SU{N) color matrices in the adjoint representation and W is the path 
ordered exponential along the x^ direction in the adjoint representation of the SU {Nc) 
group 



X ,X_l) 



Pexp 



ig / dx+A^{x'^,x ,xj_)Ta 



(2.9) 



This form is explicitly gauge invariant under the residual gauge transformations with 
gauge functions which do not depend on x~ and vanish at x^ —> ±oo. Requiring F[p\ 
to be gauge invariant, we also restore gauge invariance of the action under the gauge 
transformations which do not vanish at x^ — )■ ±oo but rather are periodic in x^ . 

This form of the interaction consistently leads to a source term in the corresponding 
Yang-Mills equation that represents classical colored particles moving along the light 
cone: 

Df.G^'" = .1+5"+ (2.10) 

with 



J+{x) = fS{x-)p\x^)t, 



TbW_ 



-[A-]TaW,+ ^^[A-] 



Expanding this expression for the current to lowest order in the field A yields 

J+(x)=.g<5(x-)p"(x^) 



(2.11) 



(2.12) 



which reproduces Eq.( |2.3| ) and is the form of the current used in As explained 

in [Q, this form is only valid in the gauge A~(x~ = 0) = 0. In more general gauges the 
current has to satisfy the covariant conservation condition 



D-J+ = 



(2.13) 



A similar Gaussian form for the statistical weight was used in in description of a large 

nucleus limit. In this case the Gaussian form is valid since the charge density is large and the color 
charges that build it up are randomly distributed in color space. 
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Our current ( ^.11 ) evidently complies with this requirement. This is a direct consequence 



of the residual gauge invariance. 

All of the above considerations finally lead us to the following effective action for the 
low X DIS 

S = i j d^x^F[p''{xi:)] (2.14) 
d'^xhj^ + J d'^x±dx- S{x-)trp{x±)W-oo.,oc.[A~]{x- , x±) 
which is the starting point of our approach. 

We end this section with a comment about the nature of the action Eq.(2.14). Although 



we use the term "effective action" when referring to it, it should be understood that it 
is different in some important aspects from the "classic" effective Lagrangians, like for 
example the chiral effective Lagrangian of pion physics. The chiral effective Lagrangian 
I/[n°] describes the dynamics of low momentum pion fields, where the momentum cutoff 
is determined by the mass of the a - particle, or alternatively the dimensionful pion 
coupling /tt, p < 47r/7r. All the modes with momenta above the cutoff, as well as all other 
heavy excitations of the fundamental theory (p, mesons, etc.) have been integrated out 
to arrive at this effective Lagrangian. In this sense our effective Lagrangian is similar. 
The gluon fields have longitudinal momenta bounded by the cutoff A+, while all higher 
momentum modes are assumed to have been integrated out. 

Importantly, the chiral physics has a sharp scale associated with it - f-^. Consequently, 
all modes of the pion field with momentum lower than the cutoff are described well 
by the chiral Lagrangian. In fact, the pions at low momenta interact very weakly, with 
the strength proportional to p'^/f-K- Therefore the perturbation theory in the chiral 
Lagrangian framework is well behaved and does not lead to large corrections to the 
tree level results. Also, the description of the low momentum pions is insensitive to the 
change of the cutoff A. 

In our case the situation is very different in this respect. There is no sharp physical 
separation scale which would separate high from low longitudinal momenta. The sepa- 
ration scale A+ we impose is arbitrary. The interaction does not die away as we go far 
below A+. Therefore there is no reason to expect that our effective Lagrangian gives an 
adequate description for the modes with momenta far below the cutoff. In fact, quite to 
the contrary as we shall see the perturbation theory in our effective theory gives larger 
corrections the farther we go below the cutoff. In this sense this effective Lagrangian is 
inadequate for description of momenta <C A+. This is of course the manifestation of 
the absence of the physical separation scalej^- This means, that if we want to describe 
low momenta, k'^ <ti A+, we have to correct the effective Lagrangian. The arguments 
presented above however fix the form of all the terms in the Lagrangian apart from F[p\. 
If the form of the Lagrangian remains the same under the fc+ evolution, the only thing 
that can change is the statistical weight e~^[''l. 

^In this respect our effective Lagrangian is more akin to "fundamental" Lagrangians of renormal- 
izeable field theories than to effective Lagrangians of the chiral physics type. 
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Physically it is quite clear what should happen. As we move to smaller values of the 
longitudinal momentum fc"*", all the gluons with momenta between and A+ are trans- 
ferred from the category of "soft" (or slow) into the category of "fast". They cease to 
be dynamical degree of freedom of interest (hence the dynamical fields have lower cutoff 
on the longitudinal momentum), but give extra contribution to the static color charge 
density p{x±). Effectively therefore as we go to lower x, the color charge density as seen 
by the soft glue changes. Since the distribution of the charge density is governed by the 
statistical weig ht e-^^P\ this means that F should change as we lower the longitudinal 
cutoff A+ . This is the physical origin of the renormalization group flow we have referred 
to earlier. In the following section we will see explicitly how this happens. First, however 
let us describe how to set up the perturbation theory in the present framework. 



3 Perturbative calculation of gluonic observables. 



The perturbation theory for the effective Lagrangian Eg. ( ^.14 ) was developed in 



It is organized in the following way. First one fixes the configuration of the color charge 
density, and performs perturbative expansion in as at fixed p{x_l). The charge density is 
not considered to be small, thus this perturbation theory is different from the standard 
one in that the calculation is performed in a non vanishing background field. In the 
second step the averaging over p{xj_) should be performed. This part of the calculation 
is contingent on the knowledge of F[p] and is completely nonperturbative. In fact the 
counterpart of this step in the standard perturbative analysis would be the specification 
of various gluon operator averages in the hadronic state. Conceptually therefore, the 
first, perturbative part of the calculation can be thought of as the calculation of the 
generalized "splitting functions" (which includes however the mixings between operators 
of different twist and is not intrinsically organized as an expansion in powers of 1/Q^), 
while the second, nonperturbative part is parallel to the calculation of operator averages 
in the hadronic state. In the standard perturbation theory, of course one does not have 
to know the operator averages in order to derive the evolution equation. As we shall 
see, the exact same thing happens in our calculation. It is only the perturbative part 
of the calculation that has to be under control in order to derive the renormalization 
group equation for F. In this section therefore we will discuss the perturbation theory 
in as at fixed p{x±). 



■ 3.1 The tree level. 

As in every perturbative calculation the first step is to find the classical soluti on to the 
equations of motion. The equations of motion that follow from the action Eg. ( 2. 14 ) are 



_9_ 



S{x-)p''{xA_)tr 



(3.1) 



As explained in the previous section, these eguations are invariant under the residual 
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X independent gauge transformation 



A + -d 

9 . 



with 



V{x) — exp [iA(a;j_. x^)] , A -^x+^±oo 



(3.2) 



(3.3) 



As a consequence, the equations of motion at fixed p have an infinite number of solu- 
tions. To properly set up perturbation theory, we should choose one of these solutions. 
Technically this is achieved by gauge fixing the residual gauge freedom. There are of 
course many possible gauge fixings. From the calculational point of view it is convenient 
to choose a gauge in which the classical solution is static {x'^ independent). It is impor- 
tant to realize that the condition of staticity of the classical solution is still insufficient. 
Even though it completely eliminates the gauge freedom of Eq. ( |3.3[ ) , there are still many 
solutions to the equations of motion. This is a consequence of the remaining gauge sym- 
metry of our problem, with gauge functions A which do not vanish at x'^ ±cxd, but 
rather are periodic in x^. With the transformation Eq.(3.3) moded out, those are 



V 



A 



-d 



9 



(3.4) 



V{x) = exp [iA(xj_)] 



To see that this is indeed the case, consider the equations Eq.(3.1) for static fields (note 
that all static solutions have vanishing A~) 



d^d+A' +g[A,,d+A'] = gp{x^_)6{x-) 
F''^ = 

The general solution to this equation has the form 

AY[p] = -e{x-)u{x^)d,u\x^) + '-e{-x-)v{x^)d,v\x^) 

9 9 
where the SU{N) matrices U and V satisfy 

a, [V^Ud.,{U^V)\ = -g^VpV^ 



(3.5) 



(3.6) 



(3.7) 



This equation obviously has a solution for any V{xj_). The matrix V, which labels these 
solutions is closely related to the gauge transformation m atrix of Eq.(3^), although this 
relation is rather subtle. Obviously any two solutions Eq.(3.6), and Af [p] are not 

related by a gauge transformation, since they solve the equation of motion with the same 
p, while the gauge transformation Eq.( [3.4| ) acts nontrivially on p. However it is easy to 
see that the set of solutions {Aflp]} with fixed p and arbitrary V is gauge equivalent 
to the set of solutions {Af [V^ pV]} with fixed V (which determines the asymptotics 
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at x~ — oo) and arbitrarily rotated p. Consequently, it would be redundant to take 
into account all static classical solutions at fixed p since we are subsequently performing 
the unconstrained functional integral over p with the measure which is invariant under 
Eq.( p.4| ). We can therefore gauge fix this extra gauge freedom by imposing, for example 
a fixed boundary condition on Ai at x~ — > — oo^. 

In this paper we will follow Ref. and choose as the subsidiary gauge condition 

d'A,{x+,xjL,x- ~oo) = 0. (3.8) 

This gauge has a nice feature that at t ^ — oo, at all finite values of z (i.e. x~ — > — oo) 
the vector potential is required to be the same as in the perturbative vacuum, Ai = 0. 
This seems very sensible, since at these times the hadron itself is still at z — > — oo and 
could not have changed the quantum state at any finite z. 



Note that this gauge condition eliminates both, time dependent gauge freedom Eq.(3.3) 



and time independent gauge freedom Eq.(3.4). It is not ghost free, and therefore the 
measure in the path integral in Eq. ( ^.14 ) must be modified by the appropriate Fadeev- 



Popov determinant 

6 {d,Aix- ^ -oo)) det (9, A [A.,{x~ -> -oo)]) (3.9) 

This modification is harmless, since as we will see later the ghosts do not contribute to 
leading order in a^, which is the order to which we calculate. 

It is important to realize that the gauge fixing Eq.(|3.S|) should be consistently used 
throughout the whole perturbative calculation. This means that not only does it deter- 
mine the classical solution we have to pick, but also the form of the propagator of the 
fluctuations of A^ around this solution to be used in the higher order perturbative cal- 
culations. In this way all potential zero modes in the propagator are eliminated and the 
calculation is unambiguous. In the standard perturbation theory, although in principle 
the situation is similar, in practice one can frequently get away without specifying the 
gauge fixing condition for the residual gauge freedom. The light cone gauge condition 
= eliminates the major part of the zero mode ambiguity and the rest of the zero 
modes start causing problems only in higher orders. It turns out that in our calculation 
we have to be much more careful, and impose the residual gauge fixing properly already 
in the lowest order. This is related to a nonstandard behavior of our fields at infinity. On 
the classical level this behavior is obvious from the form of the solution of the equations 
of motion Eq.( |3.6| ), which does not vanish at x~ ±oo. We will come back to this 
question in Section ^. 

Returning to Eq.( |3.6D we see that in this gauge for a generic fixed p{x±) there is a 
unique solution of the form 

A, = eix~)a^ixi_) (3.10) 

ai{x±) = -U{xj_)diU^{xj_) 
9 



^Alternatively, one could impose a gauge condition on p, by requiring for example that p be a 
diagonal matrix. Our choice here is dictated by calculational simplicity. 
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with the matrix U{x±) determined by 



ditti = -gp 



(3.11) 



Any gluonic observable in the tree level approximation is calculated as 



< 



0[A^] >= / D[p]e-^W O [A- =0,Ai = e{x-)ai[p]] 



(3.12) 



For example, the unintegrated gluon density defined as 



g{x,k^) =< al^{x,k^)a1{x,ki_) > (3.13) 
where a and are the light cone gluon creation and annihilation operators, is given by 




Here the < > denotes averaging over p with the weight exp{— F}. 

This has a simple reprcsciiitation in terms of the standard Feynman diagrams. The 
classical field is given by the sum of the tree diagrams for one point function in the 
background density. Using curly lines to represent gluons we have the following graphical 
representation for the full classical solution b and the first few terms in its pertrubative 
expansion 



The distribution function (up to some simple kinematical factors) is just the square of 
the field averaged over p. To the order p^ it is related in a simple way to the color charge 
density correlation function 



gp 



b = 





(3.15) 



gp 



g{x,k±) = -r^- < P{k±)p{-k±) > 




(3.16) 



where we have drawn the factors of p as diagonal lines to indicate that they are always 
associated with cikonal lines along the x~^ direction that correspond to the worldlines 
of the fast particles they represent. 
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We stress that our goal in this paper is to perform the calculation to all orders in p 
in the first order in ag. Hence an expansion in powers of p as in Eqns. ( |3l5| ) and 
( p. 16 ) would not be sufficient for our purpose. However, the above representations are 
helpful in visualizing the physical mechanism underlying the running of the charge 
density distribution with 1/x. Also, even though our interest is in the phenomenon of 
shadowing and saturation, which occur at large p, our calculational procedure should 
be valid also at small color charge density. In this limit we should recover the known 
perturbative results, which in the present context is the BFKL equation. Expansion to 
leading order in p of our result will therefore be an important consistency check in the 
calculation. 



■ 3.2 The first order perturbative corrections. 



One prominent feature of Eq.( |3.14[ ) is the full tree level x dependence of the gluon 
density. It is precisely the same as in the leading order in the standard perturbation 
theory. We know, that in the standard calculation this lowest order x dependence feeds 
back through the higher order graphs and leads to large perturbative corrections at 
small X. We expect therefore that the same will happen in our perturbation theory. 
Indeed, consider for example the graph on Figj^b, which gives one of the contributions 
to the gluon density at order a^. 





(a) 



(b) 



Fig. 1: Diagram contributing to the gluon distribution at lowest order in as (but to all orders 
in p) (a) and a typical order Us correction (b). The horizontal line represents a propagator in 
the presence of the full, p-induced background held 



This contribution was discussed in and it was shown there that it is indeed of 



order asln(l/a;) relative to the leading order result Eq.(3.14), or in the low density 



limit, Eq.(3.16). The reason for this enhancement is that when the momentum on the 
external leg l^ is much smaller than the maximal longitudinal momentum allowed in 
the field, there is huge phase space available to the emitted gluon l'^ < < A+. The 
phase space integral f then gives the logarithmic enhancement factor. 
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(a) (b) 

Fig. 2: The small p limit of the diagrams shown in Fig. ^; The lowest order terms shown in (a) 
correspond to those in Eq. ^.K^ . (b) is the small p limit of Fig^b with all gluon propagators 
perturbative. 



Carefully examining the In 1/a; corrections of Figs. |l|b, ^ , we see that it looks very 
similar to the tree level diagrams of Figs, [l^, ||a, except that the soft gluon is emitted 
not from the original charge density p as depicted in Fig.^, but rather from a modified 
charge density which in addition to p contains one extra gluon. One therefore can think 
of it as being emitted from the modified vertex of Fig.pD. Since the large correction comes 
from the region ^+ <C fc^, the emission from the modified vertex is also eikonal. So the 




(a) (b) 

Fig. 3: Original vertex (a) and vertex modihed by the emission of an additional fast gluon (b). 
"+" components of momenta are ordered from top to bottom 



source for emission of very soft gluons with Z+ <^ A+ effectively has been modified. This 
modification has to be taken into account if we are to describe properly the soft glue 
distribution. Fortunately, the change in the charge density is slow (logarithmic) so that 
this is a perfect situation for the application of the Wilson renormalization group ideas. 
We can integrate out the fiuctuations around the classical background perturbatively, 
gradually lowering the longitudinal momentum cutoff on the remaining dynamical 
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degrees of freedom. This will generate the effective Lagrangian below the new cutoff 
scale with modified F[p] . So long as we keep the change in the cutoff in every step of the 
RG small enough so that the correction to F is small relative to F itself, the perturbative 
procedure is justified. The condition for that is In ^ ^f.^ ~ 1, In ^ 1. In the 

next section we describe in detail how to set up this renormalization group procedure. 



4 The low X Wilson renormalization group. 

Let us introduce the following decomposition of the gauge field: 



Al{x)^bl{x)+6Al{x)+al{x) 



(4.1) 



where 6J^(a;) is the solution of the classical equations of motion Eq.( 3.10| ), 5A^j(a;) is the 
fluctuation field containing longitudinal momentum modes such that A+ — (5A+ = 
A"*" < < A+ while a is a soft field with momenta fc+ < A+ . Our aim is to integrate 
out the fluctuation field 6A^ in the path integral and compute the effective action for the 
soft field This integration is performed within the assumption that the fluctuations 
are small as compared to the classical fields 6°. More quantitatively, this requires that 
the coupling constant is small as 1 and at each step of the renormalization group 
procedure the ratio of the two cutoffs is not too big. In ^p- ^ ^ 

To leading order in the coupling constant we should only keep the terms up to second 
order in the fluctuation field (5A^( in the expansion of the action around the classical 
solution 6J5 (a;). 

S = -jGiaf - hA^[D-\p)rSA, + ga-J+' + 0{{a-f ) + iF[p\ (4.2) 

The inverse propagator of the fluctuation, [D^^(/9)]^'' has a nontrivial de pen dence on 
the color charge density. Its explicit form is given in the next section, Eq.( |5.2| ). 

We have introduced the modified color charge current J+', whose explicit form in terms 
of the fluctuation fields is 



J+' = p{x^)5{x-) + 5J+ + SJ+ 



(4.3) 



with 



SJ+''{x^,x+) = 6{x-) 



if'^'alSAlix- = 0) 



frV(^±) 



e{y+ -x+)-e{x+ -y+) 



(4.4) 



5A-'{y+,x- = 0) 



and 



(5J2+°(x) = -r'"'[d+5A\(x)]5A1{x) 



N, 



■p\x^)5{x' 



(4.5) 
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X / dy+SA-''{y+,xj_,x-' = 0) / dz+SA-'^{z+ ,x±,x- = 0) 



9{z+ - v+)e{y+ - x+)ivT''T''T'^T^ 

+e{x+ - z+)e{z+ - y+)trr"r''r"r'^ 

+e{z+ - x+)e{x+ - y+)trT''T'^T''T'' 



The first term in both 5 Jj^" and SJ2°' arises from the expansion of in the action while 
the rest of the terms proportional to p{x±) are coming from the expansion of the Wilson 
line term. The various terms with functions correspond to different time orderings of 
the fields along the Wilson lines. Since the longitudinal momentum of is much lower 
than of SA, we have only kept the eikonal coupling (the coupling to only), which 
gives the leading contribution in this kinematics. The contributions to SJi and dJ2 are 
depicted in Fig.^ and ^ respectively. Obviously the first diagram in Fig.^ is nothing but 
our modified vertex of Fig.pD, now cast in a more precise language. All other terms are 
nontrivial consequences of the presence of a background of fast "classical" particles that 
are encoded in the source term p and a careful treatment of the path integral over the 
modes in the interval [A+, A'+]. 




Fig. 4: Diagrammatic representation of 5Ji in terms of classical and Buctuation fields. The 
coupling to a 5y4^ field lias been indicated by a curly line whereas slow modes have been 
symbolized by dashed ones. 



We have not written out explicitly higher order in a~ terms in the effective action. There 
are of course such terms, which come from expanding the Wilson line part of the action. 
Disregarding these terms gives the effective action with the coupling of the field to 
the charge density of the form a~ J"*". However, imposing gauge invariance on the final 
result together with the requirement that the linear in a~ term of the gauge invariant 
action should coincide with the result of our calculation, the full gauge invariant form 
of the effective action will be recovered. In the following therefore we will concentrate 
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5A 




UMwmv 

SA + different time orderings 



Fig. 5: Diagrammatic representation of 5J2- Symbols as in Figl^ 



on the linear term a~ J'^ only. Note that the first term in Eq.( |4.5| ) does not have an ex- 
plicit factor of 5{x'^). However we are only interested in its low longitudinal momentum 
components since it couples directly to in the effective action. In momentum space 
this contribution is given by /°'"^/ dq^[q^ 5A\{q^)]5A^{-'q^ + k^). Since the leading 
logarithmic contributions comes from the region ^ , to this accuracy this expres- 
sion does not depend on and can be therefore approximated by S{x~) in coordinate 
space. We then can define the modified surface color charge density by 

6J^{x±,x^) — 6p{x±)d{x^) 

6p{xj_)= ldx-SJ+{xj_,x-) (4.6) 



Formally 6p defined in this way is a function of x'^ as well as x±. However, it is a 
function of 5A's which only have longitudinal momenta much larger than the momenta 
in the soft field a. The (light cone) time variation scale of Sp is therefore ^ ~ ^ and 
is much larger than the typical time variation scale of the on shell modes of the field 
a. From this point of view dp is therefore for all practical purposes (light cone) time 
independent. Technically this means that whenever we will need a correlation function 
of 6p^s, we will expand it to leading order in the time derivatives 

< 5p{x^,x+)Sp{y^,y+) >=< Sp{x^, x+)Sp{y^, x+) > +... (4.7) 

Corrections to this approximation are of order /k^ ^ x. We will therefore not indicate 
the time dependence of p explicitly. 

The procedure now is the following. We first introduce the variable p' in the path integral 

by 

i?[p,M^,a^]e^^['^-^^-''] = I D[p\p,SA^,a^]S{p' - p- dp[6A])e'^^--''^-'P^ (4.8) 



Here (5/9[(5A] is the functional of the fluctuation fields defined by Eqs.(4.4,4.5,4.6). Now 
we first have to integrate SAf^ at fixed p, and then integrate over p. 

This procedure generates the new effective action which symbolically can be written as 



exp{iS[p', a^]} = exp{-F'[p'] - ^G^ia) + igap'} 



(4.9) 
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with 



exp{-F'[p']} = J D[p, SA] 5{p' -p- 5p[5A]) exp{-F[p] - '-5AD- 



'WA] (4.10) 



Of course, to leading order in In 1/x only terms linear in In 1/x should be kept in F' . 

The integration over the fluctuation field 5Af^ is the most technically involved part of 
this procedure. We will describe in detail this part of the calculation in the next section. 
The structure of the result is however easy to understand from a simple counting of 
powers of the coupling constant as ■ Consider integration over the fluctuation field SAf^ 
at fixed p. The counting of the powers of as is done most conveniently after rescaling 
the fields and the charge density in the following way|^ 



A,, 



6p 



-A 

1 

r 

-^Sp 

r 



(4.11) 



Explicitly for the rescaled charge density we have 



p{x±_) + 5p+ + 5pl 



(4.12) 



witlT0 



5pr{x^) 



-2r^-a\5A'i{x- - 0) - -r^-p\xi_) 



X / dy"' 



e{y+ - x+) - e{x+ ~ y+) 



(4.13) 

SA-^{y+,x±,x- =0) 



and 



Sp+^ix) 



-r'^ J dx-[d+dA'^{x)]SA'r{x) 
'^p'(^±) J dy+6A~^{y+,x^,x- = 0) 
X Jdz+SA-''{z+,x±,x- ^0) 



(4.14) 



-)eix+-y+) 



(4.15) 



^The reason for this rescaling can be traced back to Eg. ( 3. 15 ). Simple counting of powers of g in 
the tree level graphs shows that for p of order all the tree level graphs are of the same order 
(see Appendix A). The classical field itself is then 0{g~^). This is also the magnitude of the field for 
which we expect to see the nontrivial shadowing and saturation effects. For parametrically smaller color 
charge densities an expansion in powers of the coupling constant automatically implies an expansion 
also in powers of p. Our primary interest is therefore in the charge densities of order aj^ . 

^We have used the identity 0{z+ -y+)e(y+ - x+) + e{x+ - z+)d{z+ -y+) + e{z+ - x+)e{x+ -y+) = 
0(2+ — to simplify the expression for (5pg. 
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with 



(4.16) 



In terms of the rescaled fields the coupling constant g disappears from the expressions 
for 6p, and appears only as the overall factor in the action. The pro pagat or of 
t he fl uctuation field is therefore of order a^. It immediately follows from Eqs.( 4.1^ ) and 
( [4.14| that 



< Sp >sA^ 0{as) 
< 6pSp >SA= 0{as) 



(4.17) 



while all other (connected) correlation functions of Sp are higher order in a^- Since we 
are working to the lowest order in as we can neglect all these other terms. Therefore 
to lowest order in a^, after integrating over SA at fixed p, we are left with the weight 
function for Sp, which generates only connected one- and two-point functions. Such 
weight is obviously a Gaussian. Introducing the following notations 



1 



< Sp^ix^) >SA-- as ln-cr''(x_L) 

X 

< Sp°'{x±,x+)Sp\yj_,x+) >sA^- asln-x''''{x±,y±) 

X 

we can write the result of the SA^^ integration in the form 



(4.18) 



D[p,p'][Det{x)]-'/^exp{~F[p]) 



X exp 



1 



2 a. In- 



Px - as In -CTa; 

X 



ly^xy ] 



Pv 



; In —c 



(4.19) 



In the above equation we adopted condensed notations: the indices x stand for the set 
of indices and coordinates {x±,a}, and repeated indices are understood to be summed 
(integrated) oveij]. 

The calculation of x and a is the subject of the following section. However, the knowledge 



of the general stru cture of the p integral, Eq.(4.20) is sufficient to perform the integral 
over p in Eq.( 4.10|) wi thout the explicit knowledge of x and a. The reason is that the 
integrand in Eq.( 4.20 ) is a function very sharply peaked around p ^ p' + 0{as), and 
the integral is calculable in the steepest descent approximation. This was done in fl^ j. 
The result is very simple 

^We note here that this result can be derived formally by introducing the variable p' with the help 
of Lagrange multiplier 

(4.20) 



5{p' -p-5p[5A]) = I D[\] 



JX{p'~p~Sp[SA]) 



and subsequently integrating out A in perturbation theory to order as 
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F' ^ F 



as ln(l/a::) 



^ S^Xuv SF SF^ 

Xuv r c c r r X^v j 

OPudpv OpuOPv Opu opy 



r^SF Sx .SF 

+Z— h z— Z- f7„ 

6pu opv opu opu 



(4.21) 



Taking the derivative with respect to In 1/a; we obtain the Wilson renormahzation group 
equation for the functional F 



dln(l/a;) 2 



Xuv 



-F 



S'Xr 



SF SF 

Xn 



dpuOPv OpuOPv OPu OPv 



,^SF Sxuv ^J(Ju ^SF 
+z- ; 1- 2- 2- — cr„ 



(4.22) 



Spu Sp-a Spu Spu 

This equation is extremely simple when written for the weight function Z = cxp{— i^} 
d 



c?ln(l/a;) 



2 dpudpv opu 



(4.23) 



Equations (4.22) and (4.23) provide the closed form of the renormalization group equa- 
tion in terms of the functionals a[p\ and x[p]- the next section we will calculate these 
two quantities. 

Eq.( [4.2£ ) can be written directly as evolution equation for the correlators of the charge 
density. Multiplying Eq.(4.23) by p^^-.-Pxr, and integrating over p yields 



dln(l/a;) 



< Pxi--Px„ > 



(4.24) 



^ ] Pxx---Pxm-lPxm + l---Pxk-lPxk+l---PxnXXmXk > 

S)<m<k<n+1 

^ < Px,---Pxi_,Pxi + ,---Px^(Jxi > 



0</<n+l 



In particular, taking n = 2 we obtain the evolution equation for the two point function 

< PxPy >= as {< Xxy + Px(Ty + Py(Tx >} (4-25) 



dln(l/a;) 



This equation is useful in making contact with standard evolution equations, since the 
correlator of the color charge density at w eak fields is directly related to the unintegrated 
gluon density in a hadron Eq. ( ji^ ) can then be straightforwardly rewritten as an 
evolution equation for the gluon density. 
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Small fluctuations in the background fleld. The cal- 
culation of a and x« 



In this section we calculate the one - point function a and two - point correlation 



function x of Sp(xj_), Eq.( 4.13|) , (4.14). First, note that these quantities are given by 
the Feynman diagrams of Fig.|6| and [7| respectively. 







Fig. 6: Diagrams contributing to a. Note that these are all virtual contributions resulting from 
contracting the SA lines of the diagrams shown in Fig.^ 



+ 



+ 



o liJiH P Miii o 



Fig. 7: Diagrams contributing to Us lnl/a;x. JVote that these arise from contracting to factors 
of the diagrams shown in FigJ^. For a separate list of all 9 contributions see Fig]^ 



The propagator lines in these diagrams are the propagators of the fluctuation fields SA 
in the non vanishing background. This is the inverse of the operator -D^J that appears 
in Eq.(U). 

At this point we see that the ghosts associated with our gauge fixing do not contribute 
to order as- The interaction of the ghost fields with the rescaled fluctuation field is 
order one, see Eq.( |3.9| ). However any insertion of a ghost vertex will lead to an extra 
fluctuation propagator and this is proportional to as- We therefore forget about ghosts 
from now on. 

Our goal therefore would seem to be the inversion of -D^J. In fact our task is a little 
simpler than that, since we only need to calculate the time independent average in a 
and the equal time correlator in x- Those are determined by the Wightman function 
of the fluctuation field Dw rather than by the Feynman propagator. The Wightman 
function satisfies the homogeneous equation 

D-^Dw = (5.1) 
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/ r 





Fig. 8: Explicit list of the diagrams contributing to x 



and is constructed from the eigenfunctions of with the zero eigenvalue. Our first 
task is therefore to find the zero eigenfunctions of D^^. For completeness we will present 
also the eigenfunctions with nonzero eigenvalues, but will not construct explicitly the 
Feynman propagator. 

For clarity we have split our calculation into three main parts. We determine the eigen- 
functions in subsection 



5.1 



find their proper normalization in subsection 5.2 and use 
these results to express the main quantities of interest, x a in subsection 
Eqs.dS^) and (|5^) 



5.3 



■ 5.1 The eigenfunctions of i^^J. 

The quadratic action for the fluctuation fields is 

5* = a^K^^ya- + 2a- {d+ Da + 2fa) + 2^+ai^-a^ + -Oi 



a, (5.2) 



Here we are using the following condensed notation 

[faf{x+,x-,x^) = r'>-S{x-)a'%x^)a''{x+,x-,x^) 
Da = A[a]a, {^^6'''' + r^''e{x-)a1)a\ 



(5.3) 



and X denotes the space time coordinates x^,a:' as well as color label. All repeated 
indices are summed (integrated) over. The function ai{x_i_) is related to p{x_\_) through 
Eq.(|]l|). The operator K is 



K 



ah 



a 
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id 



+ \2 cab 



2M''''S{x- 



where we have defined 



(5.4) 



(5.5) 



Note, that M is a color matrix locally defined in the transverse and frequency space, 
which does not depend on x^. 

For simplicity we will temporarily omit the factor l/.g^ in front of the action, remem- 
bering to restore it in the expressions for the charge density by appropriately scaling 
the small fluctuations propagator. 

We have changed our notations from the previous section and denote the fluctuation 
field by rather than (5A^. We should remember that the fluctuation fields contain 
longitudinal momenta only above some scale A"*" . The question how exactly to impose 
this cutoff is unimportant in the leading logarithmic approximation. We find convenient 
to introduce it through the infrared cutoff in coordinate space. The longitudinal coordi- 
nate x~ in our expressions therefore varies between —L and L. Whenever it is harmless, 
we will take the limit L ^ oo, which corresponds to the big cutoff ratio A+/A+ ^ 1. 

Rather than writing down the eigenvalue equations for the quadratic Lagrangian 
Eq.( ^.2| ) it is more convenient first to explicitly decouple the a~ field. This is done 
by completing the square in Eq.(5.2) 



K-\d+Da + 2fa) 



K. 



xy 



K-^{d+Da + 2fa) 



(5.6) 



d+Da + 2fa 



K- 



d+Da + 2fa 



2d aid a, — 



Defining 

a- =a- + K-\d+Da + 2 fa) 
we see that it decouples from a^. Its correlator is given by 



(5.7) 



k: 



(5.8) 

The correlator of a~ is then easily calculable once we know K^^ and the correlators of 



The calculation of if ^ is straightforward and is given in Appendix A. The result is 

. 1 , 



K-^^l - -\x- 



y 







r]{x±) \ 


\x-\ 


+ \y-\ 


2M J 



>S{xj_,yj_) 



(5.9) 



The color matrix 77°'' (xj^) projects onto the nonzero eigenvalue subspace of M. Together 
with the complementary projector fi it satisfies the relations 



liM = 0, riM = M 



(5.10) 
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and 



We note, that the operator K Eq.(5.4) has zero modes of the form 



(5.11) 



(5.12) 



and is therefore strictly speaking non invertible. The result Eq.(5.9) was obtained by 
excluding the zero mode s and inverting K on the space of functions which does not 
include the funct ions Eg. ( 5. 12 ). The normalizable zero modes oi K can not be completely 
neglected in Eq.(5^). Expanding in the basis of eigenfunctions of K 



a = I dXa^fx 



(5.13) 



we see immediately that Uq drops out from the first term in Eq.(5.2) but not from 
the second term. As a result Oq does not decouple from and the Eq.(5.7) should be 
slightly modified. In addition to the term quadratic in a" we have 



-- d+Da + 2fa 



K 



dx 



d+Da + 2fa 
d+Da + 2/a 



+ d^aid Gi - -Qi 



(5.14) 



Note that Oq does not depend on x since the zero mode of K is constant in x . The 



linear term in a~ in Eq.( ^.14 ) is in fact nothing but the Gauss' law constraint which 
remains after integrating out the a~ component of the vector potential. As we stressed 
before, our effective Lagrangian is gauge invariant under the residual x~ independent 
non-Abelian gauge transformation. As a result, the Lagrangian expanded to second order 



in the fluctuation field, Eq.(5.2) preserves the linearized version of this transformation. 
It is in fact straightforward to check that Eq.(5.2) is invariant under 



ai + Di[a]X{xj_,x^), 



+ d X{xj_,x^) 



(5.15) 



with -D^Q^] Eq.( p.3D , provided \[xj_,x^) -^x+~^±oo 0. The x^ independent part of 
a" imposes the Gauss' law constraint that corresponds to this transformation in the 
Lagrangian Eq. (E 



dx Kxytty 



{d+Da + 2/a)^ 







(5.16) 



Decoupling a 



is of course equivalent to integrating out a 
in terms of a 



from the path integral, 
except for the component of this 



This procedure solves Eq.(5.16) for a" 
equation which is proportional to the zero mode of K , since this component does not 
contain a~ . This component of the equation is a constraint that involves Oj only and 



should be kept intact in the path integral for a^, Eq.(5.14). The field Oq is just the 
Lagrange multiplier that imposes this constraint. 
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Now that we have disposed of a~ we have to find eigenfunctions and eigenvalues of the 
operator D~j^[p\ defined by the action Eq.(5.14). It is convenient to parameterize the 
fields ai in the following way 



+ e{-x-) 



a''_{xj_, ,x^ ,x )+Y-{xa_,x^) 



(5.17) 

The reason we choose to use this parameterization is that the equations of motion 
(the eigenvalue equations) as derived from Eq.(5.14) are first order in and contain 
coefficients of the form d{x~ ). We therefore expect the eigenfunctions to be discontinuous 
at x^ — 0. Also, since the classical background fields do not vanish at x^ oo, we 
should allow the same asymptotic behavior in the fluctuations. We have separated out 
for convenience the components of the field jj. which do not vanish as x~ ±oo so 
that by definition 







(5.18) 



Substituting ( ^.17 ) into the action ( |3.14 ) we obtain 



S 



dx 



d-a\d+a\ + -a\Dla] 



dx 



d-ald+c 



1 



[7; + v\_]d-[Y_ + 1-1] + [d-f_v'_ - d-Y+v\_] 



1l 

2 



d-' 7:^ — Z?-* 7^ — a' [vi^ + 



where 



= a±(a;" = 0) 

The covariant derivative in this equation is 

Df = {8^5"'' + r'"'a^,) 



(5.19) 
(5.20) 



(5.21) 



We hope that the use of the same symbol as in Eq.(5^) does not cause confusion 
In this parameterization the linearized gauge transformation acts as 
dY+ = D'A, S-fL = d'A, Sa^ = d-A 



(5.22) 



The equations for eigenfunctions are 
SS 



Sai+ 



0ix-) 



2d+d- 



a\_ + S{x-) 



9- [71-7^ 



+a'— a [dj+ - D-f- - a(v+ + V-)] + a"^ fia^ 
oa 



(5.23) 
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SS 



9{-x- 



2d+d- + dl 



a'_ + S{x ) 



(5.24) 



SS 



oa 



(5.25) 



SS 

S^i- 



oa 



= AL7! 



(5.26) 



where all the derivatives are with respect to transverse coordinates unless explicitly 
specified. These equations are supplemented by the constraint 



(5.27) 



First consider the zero eigenvalue A = 0. Due to the gauge symmetry, the equations for 
eigenfunctions have infinity of solutions for A = 0. However, as stressed in Section |^ wc 
must work in a completely fixed gauge, which we have chosen as diai{x^ —00) = 
0. In the notations of this section this means d'j- = 0. With this gauge fixing it is 
straightforward to find the solution 



d^Pt 



-X ) exp ( - ^Vtxt ) wl,r(Pt) 



+ 6'(x-)C/^(xj^)exp (^i^x- - iptxt^ [Uv+^r] (P-l) 
+ 0(-x-)7l,.(x^) + 0(x-)7;,. 



(5.28) 



The frequency is a good quantum number since our background field is static. Here r is 
the degeneracy label, which labels independent solutions with the eigenvalue A = and 
frequency p~ . In the free case it is conventionally chosen as the transverse momentum, 
{r} = {pi}- The matrix U{x±) is the same SU(N) matrix which defines the classical 
field Eq.(3.10) The auxiliary functions 75-, f± are all determined in terms of one vector 
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function. We take this independent function as wLH. Then 



Si] _ 2D'—D^ 

D2 



(5.29) 



For the eigenfunctions corresponding to eigenvalues A 7^ there is no gauge invariance. 
Accordingly the functions vanish at infinity j±,\^o — and the solutions are 



<fpt 



9{-x ) exp ( i X - iptxt ] v'^^x.riPt) 



+ 6{x )C/^(a;_L)exp '^^^ ^* x - iptx^ {Uv\ .^ x{Pi^)) 



(5.30) 



with 



a — a(w+ + u_ 



(5.31) 



We now have to construct a complete set of solutions, which is tantamount to picking 
rC^^J-) foi' every eigenvalue A as a complete basis of functions on the plane. This basis 
should be chosen such that the solutions Eq.(5.28) are properly normalized and are 
orthogonal for different values of r's. 



■ 5.2 The normalization of the eigenfunctions. 

The orthonormality condition for the eigenfunctions a\ p_ ^ i^ 

j d^xaX^^^,ix)a\\^^,^A^) = S{X- X')S{p- -p-')S{r-r') (5.32) 

Although for the purpose of our calculation we only need eigenfunctions with the eigen- 
value A = 0, it is convenient to consider the orthonormality relation for arbitrary A. The 
reason is that if we take A = A', the factor 6{X — A') gives a divergent constant, and it 
is difficult to determine the numerical coefficient in front of it. Taking A and A' generic, 
we can explicitly extract the S - function factor and determine the coefficient. 

* These expressions are valid up to terms of order 1/L. The omitted terms do not contribute to the 
leading order in Inl/x. 

^One could ask whether the presence of the Lagrange multiplier in the Lagrangian can modify 
the normalization condition for the eigenfunctions. It is sh own i n Appendix B that this is not the case, 
and the appropriate normalization condition is indeed Eq.(|5.3^). 
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Let us consider the scalar product 



2p- 



where 



\- \' +i€ 

+ L{N_ + N+) 



2p- 



X- X' -ie 



-M4 



(5.33) 



M+ = J d?x^v^^^^,{x^)v^:,,,^,{x^) 

N- = j d'x^Y-\A^^h"':,,^y{x^) 
N+ = J d'x^^'^+\,{xi_hr:r',x'i^^) 

here again we have kept the terms of order and dropped the cross terms since 

dx'a'Y l/L (5.35) 



(5.34) 



and this can be ignored at large L. 

Consider first the case when both eigenvalues A and A' are non-zero. It follows from 
Eq.(|3|) 



where S is the operator 



da 



This operator satisfies, 

[S, S^] ^ 
and therefore S^^S^ is unitary, so that 

V^V-\- — V_V- 



(5.36) 
(5.37) 
(5.38) 
(5.39) 



and A/+ = M_. The orthonorniality condition ( 5.32 ) then becomes 

d^'x a\-^^-^^^{x)a\\j,-,yix) = 27rS{X - X')6{p- ~ p-')2\p-\M_{r,r') (5.40) 



It is clear now that for A 7^ 0, we can take our orthonormalized basis to be 



1 



\J\Tip 



A^^O 



(5.41) 
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For A = 0, there will also be a non-vanishing contribution from the term which involves 



7_i_. The x~ integral in the normalization condition Eq.(5.32) gives a factor of L. Com- 
paring it with Eq.(5.4C) at A = A' we identify this factor as 47r|p~|(5(A — A'). We are 
then left with 



M_ + M+ + N^+N+ = 



1 



47r|p 



j6{r - r') 



(5.42) 



It is easy to see that for the zero modes the relation between u_ and w+ is also unitary. 
Using M_ = M_|_ and the explicit expressions for 7+ we get 



Att\p 



5{r-r')= / d^x^d^y^v%{x^)0'^{x^,y^)vl^^,{y^) 



(5.43) 



with 

0*^'(a;±,2/i) 



ski 



yz JJZ 



jjZ 



(a;±,y±) 



(5.44) 



This relation means that for proper normalization one should choose the functions w!_ ^ 
as eigenf unctions of the operator O The normalization of v'[_ ^ should not be one, but 
rather \/Vk where k is the appropriate eigenvalue of the operator O. The degeneracy 
label r therefore numbers the vectors of this particular basis. Since O is a Hermitian 
operator, its eigenfunctions form a complete basis, and therefore the basis of our eigen- 
functions is also complete. Therefore we have 



4:Tr\p 



All of our results then will be expressed in terms of O ^ where 



1 + 2[(1 - d^D)il - D^d) - (1 - D^D)] 



-Ikl r 



(5.45) 



(5.46) 



Q2 



■ 5.3 The induced charge density. 

We are now ready to calculate the induced charge density 5p. As was mentioned in 
the beginning of this section, since we are interested in the equal time correlations of 
the fluctuation fields ai we will only need the on shell propagators, and therefore only 
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eigenfunctions for the eigenvalue A = 0. To see this exphcitly let us consider a typical 
expression we have to evaluate in order to calculate the charge density correlator 



dx~ 



dX 



0)> 



X + ie 



dp d k±d p± e 



.2X + p^, +k 

I i 



2p- 



(5.47) 

' F(k^,p^)vl^,^^{p^)vl*{k^)d{x~) 



We have regulated the A dependence of the integrand by moving slightly away from 
x^ — 0. At nonzero x~ we can close the integration contour in the A plane. At every 
fixed value of p~ the contour can be closed either in the upper or in the lower halfplane, 
depending on the sign of x~ . The only contribution to the integral comes from the pole 
at A = 0. If the contour is closed upstairs the integral vanishes, while if it is closed 
downstairs there is a contribution 27ri. Therefore for every p~ the A integral gives a 
factor 2'Ki9{x~) or 2iTi9{—x~), depending on the sign of p~ . Multiplying by 5{x~) and 
integrating over x~ gives the factor 1/2 in either case. The result of the A integral is 
therefore that it puts the propagator on shell (A = 0) and gives the numerical factor ni. 

In the following formulae, A = is therefore assumed. 

First, we calculate . In fact as is obvious from the explicit expressions for the charge 
density, Eqs.(4.13) and (4.14) we need only a^{x^ — 0). Also, as can be easily checked 
oT d oes n ot contribute to the order In 1/a;, and we omit it in the following. Then, using 



Eqs.dl^ll), (1^3) we find 



d+Da + 2fa 



ip-e'P ^ 



Z?2 



{y- 



(5.48) 



r(a;_L) 



Q2 



fp ,r 



the integration over p^ and r is implied in this equation. The objects fp-,r are the 
coefficients in the expansion of the fields in the basis of the eigenfunctions of the 
operator 

aii^) = JdXdp~dr ai^x,p- ,r{x)f\,p- ,r (5.49) 
Since our eigenfunctions are properly normalized, /'s have standard correlator 



< fx,p -rfy p-' 



X + ie 



-5{X- X')6ip- -p-')Sir-r') 



(5.50) 



Now, using the Eqs.( ^1.13D ,( |4.14| ) and ( |5.2S| ) we find 



fp ,r 

(5.51) 
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and 



Dv+ 



Dv* , dv* , 



5? 



(5.52) 



f 



abc 



/O /"oo 
-oo ' ' ' ' Jo 



+ ,p-' ,r' 



P ,rJp-',r' 

Now it is straightforward to evaluate a and x- Since Spi is linear in Ui, clearly < Spi > = 
0. Also, 6pi is order g while Sp2 is order . Therefore to order as only Spi contributes 
to Xj and only 6p2 contributes to a. The frequency integral is trivial. The only p~ 
dependence is in the normalization factor Eq.( ^.45" ). The integral over p~ then 

gives the logarithmic factor which we identify with Inl/x. 

Using the explicit expressions for w+ and 7+ in terms of w_ and the normalization 
condition Eq.(5.45) we obtain 



(5.53) 
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r)i fik ■ 



(x±,z±) 



kn 



gh 



X < a 



In 



fh 
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de 
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I fj.mn imn\ 



fe 



jacejbdepb^^^^ 
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n eg 



{x±,yi_) 



7-)fe al 



dh 



{x^,yi_) 



(5.54) 
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where R^'^{x±) is 



<PP1. A± „2„-ifa,-.,)x, 



(27r)2 (27r)2 

and we have defined the projection operators 



(5.55) 



QIQ3 



(5.56) 



These expressions can be somewhat simplified. The inversion of O can be performed 
exphcitly as far as the transverse index structure is concerned. 



[(t-l)0-\t-mnx^,y^) 



(5.57) 



< xi_\5: 



ah 



ah 



Here < > denotes the configuration space matrix element in a standard way. 

We also find it simpler to use the matrix notation 



The rotational scalar operator S is defined as 

1 ^,Di di .,Di di . 
O — I of - -\( . 1^ 



(5.58) 



(5.59) 



In terms of this operator we have 



X'"'{x^,y^) = -<x^\ 



cr'^ix^) 



g. £)■ 
2ai-{a,d}-^-{a,D}^ 



yQ2 2)2- 



^Q2 D2' 



2aj-^{a,d}-^{a,D} 



ah 



■ a D, _ D, 



\y± > 



(5.60) 
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Here {X, Y} denotes anticommutator. 
Equations ( ^.54 ), ( 5.53 ) and ( 5.6l| ) are the central result of this paper. Those are expres- 



sions for the coefficient functions of the renormalization group equations as functions 
of the charge density p. These expressions do not look very simple and clearly one will 
need to develop some intuition and deeper understanding to be able to use them in 
either analytic or numeric calculations. This is the matter of future work. For now we 
are content with being able to derive these explicit expressions. 

As an important cross check on our results, we have checked that in the weak field limit, 



expanding the renormalization group equation Eq.(4.22) to leading order in the charge 
density p we recover the BFKL equation. The details of this calculation are given in 
Appendix C. 



6 Discussion 

The main result of this paper is the calculation of the one - and two - point functions 
of the charge density induced by the gluonic fluctuations. This completes in the formal 
sense the derivation of the renormalization group equation that describes the flow of 
gluonic observables at low x according to the ideology of [0-|l3|. We want to point 
out that in fact, the flow is described not by one RG equation, as in a theory with one 
running coupling constant and not even by a flnite set of equations, as in a theory with 



finite number of relevant operators, but rather by a functional equation Eq.( 4.22 ). The 
functional equation is equivalent, of course to an infinite number of ordinary equations. 
This can be interpreted as indicating that the low x RG flow has an inflnite number of 
relevant operators. This is a rare example of the renormalization group flow in an inflnite 
dimensional space of relevant couplings with all "/3 - functions" calculable explicitly. 

Much work remains to be done to understand the physics of the full nonlinear evolution 
equation. It is probably wise flrst to see whether one recovers the simpler known equa- 
tions as its particular limits. As we have mentioned, we have checked explicitly (see [ [l^ 
and Appendix C) that in the leading order expansion in powers of the charge density 
our equation reduces to BFKL equation. The DLA limit of the DGLAP evolution is 
also obtained if we expand to leading order in p, impose the transverse momentum or- 
dering in the rungs of the ladder in the real diagrams of Fig.^] and neglect the virtual 
contributions of Fig.^. With a little more work one should be able to recover the GLR 
equation Q. To this end one has to expand our result to the next to leading order 
in the charge density p and impose the DLA kinematics. Without imposing the DLA 
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transverse momentum ordering the next to leading order expansion should reproduce 
the triple pomeron vertex. 

These are important consistency checks on our calculation and they should certainly be 
performed. In the framework of the full nonlinear problem, there are two very interesting 
questions which can be asked immediately. First, does the flow described by Eqs.( 12^ ), 



( p. 54 ) and (5.53) have a fixed point. The presence of such a fixed point would mean in 



our framework the saturation of gluonic observables at low x. Needless to say, if such 
a fixed point exists and the fixed point value of F[p\ can be determined, it would be 
extremely interesting. It would describe a universal behavior of DIS observables at low x, 
independent of the hadron that is being considered. The statistical weight e'^^^^ defines 
a two dimensional Euclidean field theory. It is interesting to note, that the evolution 
equation itself does not contain any scale. Therefore, if the fixed point exist it could 
be scale invariant, and in that case very likely also conformally invariant^. It may 
be therefore possible to study it with the methods of two dimensional conformal field 
theory. 

Second, even if the fixed point does not exist it would be interesting to investigate what 
is the impact of the low x evolution on observables at small transverse momentum. 
Starting from some reasonable initial condition at xq and evolving F to low enough 
values of x, one could study the low k± behavior of the resulting two dimensional theory. 
Again, it could be that at low k± the two dimensional theory becomes conformal and 
can be analyzed analytically. 

Another outstanding question, is how to generalize this approach to include not just 
DIS but also hadron - hadron collisions. An effective Lagrangian for the hadron collision 
in multi Regge kinematics was derived by Lipatov Q . It would be worthwhile to extend 
the renormalization group approach to this case. We note that the leading order of the 
perturbative calculation for two hadron collision was considered in p7[ | to first order in 
p and numerical work is in progress to understand the nonlinear effects to leading order 
in a, @. 

We want to conclude with a discussion of the physical picture of the difference between 
the BFKL limit and the low x DLA DGLAP limit as it emerges from our approach. 
Although this does not have direct relation to the nonlinear problem considered in 
this paper, the observation can hopefully help to put our approach in a more general 
perspective. 

It is interesting to interpret our calculational procedure from the perspective of the Born 
- Oppenheimer approximation. The Born - Oppenheimer approximation is standard in 
systems that have two distinct time scales. One considers the slow degree of freedom X 
as a static background and solves the dynamics of the fast degree of freedom Y with 
given background X. Integrating out Y generates a change in the Lagrangian for X. 
This is of course the standard procedure for deriving effective Lagrangians in theories 
which contain well separated fast and slow degrees of freedom. An example is the chiral 
Lagrangian, where light pions are slow, and heavy p, (j) etc. mesons are fast. From this 
point of view, in our case we would like to think of the partons with large longitudinal 

^"We thank M. Wiisthoff for this observation. 
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momentum p'^ > A"*" as the slow modes, since their frequency P~ ^ is small|^ The 
charge density p due to these partons is considered to be static while integrating out the 
fluctuation fields SA with momenta A+ > A:+ > A + . However as stressed before, our 
system does not have two sharply separated time scales, but rather a continuum of time 
scales. We are integrating over SA in order to derive the Lagrangian for the soft fields 
a. Those contain light cone frequencies l~ ^ ^ ^ . The fluctuation fields 6 A would 
appear as fast modes relative to p but as slow relative to soft fields a. The situation in 
fact is slightly more complicated, since in principle the field modes SA contain also vastly 
different transverse momenta. For k'^ > p\ they are in fact faster than p, however for 
very small transverse momenta k^ <C p\ the frequency of the 5 A field can be as small 
or even smaller than p~ . The splitting in terms of longitudinal momenta fc"*" therefore 
does not exactly correspond to the Born - Oppenheimer type splitting in terms of th?e 
frequency. 

The charge density Sp that is induced by the fluctuations SA accordingly contains two 
types of contributions. First, there are contributions with frequencies of order k~ 3> 
p~ . Those come from the real diagrams of Fig.^ with ordered transverse momentum, 
^ Pj_- For convenience of reference we redraw one diagram from this set in Fig.0 

The frequency of this component of the induced charge density Sp faster is ^ 

There is another component Spsiower which comes from the unordered region of the 
transverse momenta fc^ ~ p^. This component is in the same frequency range as the 
original p. Since it is not contributed by faster modes it is not a Born - Oppenheimer 
type contribution. In principle there is also a contribution from the region with oppo- 




Fig. 9: Part of the real contribution to the induced color charge 

site ordering /c^ ^ p^. This component of the induced charge density would contain 
frequencies which are even smaller than those in p. However it is exactly canceled by 
the virtual corrections of Fig.^ as can be checked explicitly from the expression for the 

^^Note that those are the modes that we called "fast" in Section ^ We hope this does not cause 
confusion. Indeed these modes are "fast" if we consider their variation in time t, since their energy is 
large. However, in light cone time x"*" = z + t these same modes are almost static, since their light cone 
time dependence is given by exp{ip~x~^} and p~ is small. In this section we are interested in the light 
cone time variation. 
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BFKL kernelQ. 

So, to reiterate, as we move to lower x we get two distinct contributions. First, there 
is a Born - Oppenheimer type contribution. It appears because we "renormalize" the 
notion of staticity by including new (faster) modes in the category of static. This is the 
faster- Physically, this relates to the processes at smaller x that happen at ever faster 
time scales. Second, there is a new contribution to the charge density which has a long 
wavelength but also low frequency. This is 5 p slower- This is not a contribution of the 
Born - Oppenheimer type and appears since the splitting in terms of the longitudinal 
momentum does not strictly coincide with the splitting in terms of frequency. 

Now, the DLA approximation to DGLAP assumes transverse momentum ordering and 
therefore includes only 5p faster in the induced charge density. The BFKL evolution 
on the contrary includes both contributions. Therefore the DLA DGLAP includes in 
the induced density only contributions which are faster than previously present, and 
is evolution simultaneously in the longitudinal momentum and frequency. The BFKL 
evolution includes contributions which at every step in the evolution modify also the 
slow component of the charge density, and is therefore evolution only in the longitudinal 
momentum. 

The Spsiower contribution to the BFKL kernel is known to be problematic, since it cre- 
ates a channel through which the nonperturbative small transverse momentum modes 
couple to the evolution. The result is the infamous random walk in the transverse mo- 
mentum space in the asymptotic solution of the BFKL equation. Our full nonlinear 
procedure is similar to the BFKL approach in the sense that it docs not discriminate 
between the field modes on the basis of their frequency. Whether it still suffers from 
the same low transverse momentum problem is not clear to us at this point. It is pos- 
sible that the nonlinearities in the equation suppress the low transverse momentum 
contributions by generating a dynamical scale somewhat like what happens in the fi- 
nite temperature and finite density equilibrium systems. It would be very interesting 
to develop a renormalization group procedure in which the evolution parameter is not 
the longitudinal momentum (like in our approach) and not the transverse momentum 
(like in the DGLAP evolution) but rather directly the time resolution scale. This kind 
of approach would consider the contributions of low k± modes at low x as part of the 
initial condition rather than part of the evolution and would thereby provide a cleaner 
separation of perturbative and nonperturbative effects. A nonlinear evolution obtained 
in this type of approach should be closely related to the nonlinear evolution equation of 
Laenen and Levin |po|| . 



^^One can be tempted to think of these virtual corrections as the Born - Oppenheimer type back- 
reaction, which "renormalizes" the effective distribution of the slow component of the charge density. 
This is not the case since the contribution in the virtual diagrams comes mostly from the low transverse 
momentum region in the integral, and is therefore due to slow modes in the loop. 
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A High and low density situations: the parametric 
size of p 

The situation that interests us in this paper is the one when the color charge density is 
high enough so that the nonlinearities are important. Aheady when looking at the tree 
graph expansion of a classical field generated by a color charge density p, we are in a 
position to judge how big p has to be parametrically for all of the terms to contribute 
equally to the result. 

To do this wc observe that one may build all diagrams that make up the classical field 
starting from the linear term 

«iij)iji!s.(*>m' gp (-^-l) 
by stripping off a factor gp and replacing it successively by either 

gp gp 



g W or g"^ iiuiiffloiuffluujj) gp 



gp gp 

To generate diagrams which are all of the same order we need to simultaneously satisfy 
the equations 

g ■ order {p) — g {g ■ order{p))^ (A. 2) 

g ■ order(/9) = 9^ id ■ order(p))^ (A. 3) 

This fixes 1 = g"^ ■ order(p) i.e. p to be of order g~'^ . 

It is therefore for p ~ 0{g~^) that the nonlinear effects described by the RG evolution 
derived in this paper are physically important. 




B Inverting K. 



In this appendix we invert the operator K which appears in the small fluctuation action 



K" 

xy 



(B.l) 



First, note that the frequency and the transverse coordinate dependence of K is trivial, 
and therefore p~ and x± are conserved quantum numbers in this inversion. For the 
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purpose of this calculation we can imagine that the color matrix M has been diagonalized 
at every point in x±^ and we will therefore treat it formally as a number. 

Let {x) be the set of eigenfunctions of K 



dy~K'^\x,y)f\p-,y-,x^) = Xr{p-,x-,x^) 



(B.2) 



Assuming for the moment that the eigenvalues are positive, the eigenfunctions have the 
form 



f^.(x)=e{-x-)e^P '^"U'^ix^,p-)+0{x-) 



(B.3) 

where F and G are to be determined. Requiring that /° is continuous at a; =0 gives 
jja _ pa _|_ Qa rpj^^ ^^^^ derivative of /" shoul d be discontinuous such that when 
acting on it will cancel the S{x~) in the equation (BJ_). In other words, the discontinuity 
in the first derivative must be equal to the integral of the delta function across the 
discontinuity. This gives 



fcont = Ke^^^ ^e{x-f^^^^^r''^U'^p^x^). 



p+p 



(B.4) 



These eigenfunctions correspond to a positive eigenvalue 

To determine the proper normalization of the "polarization vector" it is easier to 
work with the symmetric and anti-symmetric combinations 



1 



fLym = U^cosp+x- - f^o^U^p^ sin p+ X 



e{x^) - 9{-x-) 



and 



J Ota 



iU^ sinp+a; 



2p+p 



The antisymmetric functions vanish at a: =0 and therefore are the same as for M = 0. 
Their normalized form is just 

faasym ~ sin(p'^a; 
The overlap matrix for the symmetric eigenfunctions is 







t 




j dx 






/r/3,p+')(^) 



5{p+-p+') 



2 1 



{2p+p-y 



and we get 



(p+)2 + M2 
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These are the eigenfunctions for "continuum" states, the ones corresponding to positive 
eigenvalues. 

Note that since trf ■ p — 0, the matrix AI has both negative and positive eigenvalues. 
Therefore the eigenvalue equation (B.2) must also have negative eigenvalues and corre- 
sponding "bound state" solutions. The bound state wave function must be symmetric 
under —x~ and has the same form as the continuum solution except that p+ is 

imaginary. Requiring that it decay exponentially at large x~ , we have 



p+ = iMB{M) 



The eigenfunction then is 



0{x-)e 



-Mx 



+ e{-x-)e 



Mx 



(B.5) 



where iV is a normalization factor and Vj^ are the set of eigenfunctions of M with positive 
eigenvalue. The normalization factor N is easily calculated and is given by N"^ = M . 

It is easily checked that the set of our eigenfunctions is complete. The completeness 
relation is 



/ = 

This can be written as 



a'^ix^p) 



a''{y,p) 



Six - y)S'^' 



(B.6) 
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I = S{xj_-yA_) I -^h{x ,y ,p 



where for the continuum solutions we have 
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+ {x^,y 

while for the bound state 
h = \i\M\+M) 

+ {x'^,y^^-3 



i(| M I -M)e{-x- + y-)el^'^l(^"-2'") + 6{x- - y-) 
-x~,-y') 



(B.7) 



e{x-)9{y-)e-\^\^-^ +y ) +6l(x-)6'(-y-)e-l*^l(^ ^ 

-.-y-) 



B. 
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All M dependent terms cancel out between the two contributions. Performing the p+ 
integration and adding all the terms establishes the completeness relation ( |B.6D . 

Now that we have the complete set of eigenfunctions, we are ready to invert the operator 
K. But first we should understand the zero modes p+ = 0. These eigenfunctions do not 
vanish at ±cx3 and we should be careful with them, for example when integrating 

by parts. In fact when calculating the propagator, these eigenfunctions should be ex- 
cluded entirely from the sum as explained in Section |5[ The convenient way to do that 
is to regulate the factor 1/A which enters in the calculation of the propagator as 



1 



A 



A A2 + £4 

taking the limit e ^ at the end of the calculation. We use the same regulator to 
regulate a possible singularity M — > in the bound state eigenvalue. 



1 



AP M4 + £4 
The propagator is calculated as 

(B.9) 

The result is 
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M 
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J\/ie(x +y ) +y ) 



M + iVie 



M 



'iy/itix -y ) p-Vif-(x -y ) 



M + iVii M-Vie 



-y ) \5{p -p )S{xi_ - vi_) 



(B.IO) 



Expanding the above expression in powers of e to order one, we get 



1, 



1 



7^ 



\y 



■q 



2M 2V2t 



5{p-~p~ )8{x^-y^) (B.ll) 



where we have defined the projection operators t] and /i that project on nonzero and 
zero eigenvalue subspaces of M respectively 



(B.12) 
(B.13) 



Note that the last term in Eq.( p3.11 ) diverges in the limit e ^ 0. However examining 
carefully equations in Section]^, we see that always acts on a particular combination 
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of fields B which satisfies the constraint fxB = by virtue of the Gauss' law. We can 
therefore omit this term from the expression for altogether, which is what we did 
in the text Eq.(O). 



C Proper normalization of eigenfunctions. 

In this appendix we show how to properly normalize eigenfunctions in a theory with a 
Lagrange multiplier field. Consider a quadratic form 

S = a^M.jaj, i = 1, N (C.l) 

where the variables Ui are constrained by m linear conditions 

c>i = 0, a = l,...,m (C.2) 



Our problem is to invert M on the subspace whose vectors satisfy Eq.(G.2). This is 
the precise analog of the system we deal with in Section ^. Let us assume that all m 
vectors cf are linearly independent. In that case they span an to - dimensional subspace 
C of the original N dimensional vector space V. Let P be an orthonormal basis on this 
subspace. We can then construct the projection operator 

R, = (C.3) 

which projects on C. 

Then instead of considering the matrix M we should consider 

M = (1 - P)M{1 - P) (C.4) 
and invert it on y — C . The eigenvalue and eigenfunction equations for this problem are 

Ma = A(l - P)a (C.5) 

or alternatively 

(1 - P)Ma = Xa 
c"a 

Clearly the eigenfunctions have to be normalized in the standard way 

a'^af = 5a,a' (C.6) 
The inverse of M on V — C is constructed as 

M-'=Y.-^ (C.7) 
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since 



MM~^ = (1 - P)MM-^ = a^a*/ = Py (C.8) 

A 

An alternative to explicit construction of the projection operator P is introduction of 
the Lagrange multipher field, as we did in Section ||. We add to the Lagrangian the term 

6"c"a (C.9) 

The eigenvalue equations we have to solve now are 

Ma + b°'c°' = \a 
c^a = a 

Now since, Pc" — c" and Pa = 0, this equation gives 

PMa + 6"c" = (C.IO) 



Solving this for 6" and substituting back into Eq.(C.lC) we obtain again Eq.( p.6D . This 
proves that the eigenfunctions obtained through introduction of the Lagrange multiplier 
are the same as in the straightforward calculation in which the constraint is solved 
explicitly through construction of P. It is clear th en th at the normalization of these 
functions should be the standard normalization Eq.( |C.6[ ). 



D The BFKL limit of the general evolution equation. 



In this section we will show in some detail how our general expressions for x"''' and ti" 
give the BF KL ke rnel. To do so, we have to take the limit of small p in the evolution 
equation Eq.(4.22). In fact it is more convenient to consider directly the equation for the 
density correlation function Eq.(4.25). As was shown in using Eqs.(3.14) and (4.16) 
the unintegrated gluon density 4'{kj_) (which is the quantity which evolves according to 
the BFKL equation) to leading order in p is just the charge density two point function 
< p{k±)p{—k±) >. The BFKL equation should therefore be just the weak field limit of 
Eq.(P^). 

To verify this we will need the expressions for x"'' and cr" expended to first order 
in p. Let us consider the contributions of the real diagrams given by Spreai in equa- 
tion (refeqitotalrhoreal). Using the expressions for v'\, and and expanded to 
first order in a 



1 



2^ 

dl 



[da 



W - 7^(c'a + «a)-2- 



(D.l) 
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leads to 



spnx±) = -2/"" 



In the momentum space representation 



9wl 



(27r)2 



d'y±- 



With these expressions we have 
To this order the normahzation of t;,- is 



P± 



v-iy±) 



dl 



We then have 



2,71 z 
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(D.2) 



(D.3) 



(D.4) 



(D.5) 



< 5p1{x^)6p\{y^) >= g'r'' r'-\nllx\a'{xi)a''\yi_)5\x^ - y^ ) 

TT 



a^%x^)p%y^)-p%x^)a^%y^) 
(Ppi_ 1 



Jp±ix±-y±) 



+ p^x^)p'^{y^) 



(27r)2pi 

This expression coincides with Eqns. (53,54,55) in |12|. Using 

Q,i = _^ we obtain 



(D.6) 



SpUk^)6pU-k^) 



(2 



^ J Sp^p''{pA_)p%-p^] 



kl 



p]_{pi_ - 



(D.7) 



This is precisely the real part of the BFKL kernel. 

It is straightforward to repeat the above procedure for the contribution of virtual dia- 
grams < 5p2{x±) > (5.52). In this case, the term i?"'' vanishes to order a. Using our 
expanded expressions for a+ and 74- and noticing that 7+ starts at order a, the first 
line in (5.52) gives 



<5p,{x^)>^,)—-f f Jdy^piy^)J^ — 



d^p±_ d^q±_ PA_-qi_ ^ 



|2 p^ql 



ip±+q±){x±-y±) 



(D., 



This agrees with the corresponding term in [g_2|. In the second line in the expression 
for Sp2 (5.52), we can take the order inside the brackets since there is already an 
explicit factor of p present. This gives 



< Sp^^ix^) >(2)= -0^\nl/xp'^{x^) 



Pi 



(D.9) 
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which is exactly equation (50) in Collecting all the contributions, substituting them 
into Eq.(4.25) and identifying the density correlator with the unintegrated gluon density 
6 we obtain 



dlni 



(27r)3 



(D.IO) 



This is precisely the BFKL equation 
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